MULTIPLICATION OF FREE RANDOM VARIABLES 
AND THE ^-TRANSFORM: THE CASE OF 
VANISHING MEAN 

N. RAJ RAO (*) AND ROLAND SPEICHER (t) 

Abstract. This note extends Voiculescu's S'-transform based an- 
alytical machinery for free multiplicative convolution to the case 
where the mean of the probability measures vanishes. We show 
that with the right interpretation of the S'-transform in the case of 
vanishing mean, the usual formula makes perfectly good sense. 



1. Introduction 

Multiplicative free convolution IE was introduced by Voiculescu 0, 
as an operation on probability measures to describe the multiplica- 
tion of free random variables. Since various classes of random matrices 
become asymptotically free, the multiplicative free convolution is an 
important concept to deal with the asymptotic eigenvalue distribution 
of products of random matrices [HI H] • A powerful analytical tool for 
an effective calculation of Kl is Voiculescu's S'-transform. In the usual 
presentations of the S'-transform it looks as if its definition breaks down 
in the case where the considered random variable has mean zero. Since 
interesting cases, like the multiplicative free convolution of a semicircle 
with a free Poisson distribution (corresponding in the random matrix 
language to the product of a Gaussian random matrix with an inde- 
pendent Wishart matrix), fall into this class one gets the impression 
that the S-transform machinery has a crucial flaw here. However, as 
we want to point out, this is only ostensive, and with the right inter- 
pretation of the S-transform in the case of vanishing mean, the usual 
formula makes perfectly good sense. We will show this by combinato- 
rial arguments in the setting where all moments exists. It stands to 
reason that one should be able to establish this using purely analyt- 
ical tools, for more general distributions for which the combinatorial 
moment based arguments do not suffice. 
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2. Main result: one mean is zero 

Let us first recall the definition of the S'-transform. We start by 
working on an algebraic level, further below we will address the question 
of positivity (which is crucial for the definition of Kl). So a random 
variable x is an element in some unital algebra A, which is equipped 
with a linear functional ip such that = 1. Our main interest is 
in the moments (p{x^) of the random variable x. Of course, of main 
interest for us are real-valued random variables X : M in the 

genuine sense, where is given by taking the expectation with respect 
to the underlying probability measure, 

¥;(X") = [ X{uYdP{uj). 
Jq 

In our algebraic frame we have to restrict to situations where moments 
determine the distribution uniquely; let us, however, remark that by 
analytic tools much of the theory of multiplicative free convolution can 
be extended to classes of probability measures with no assumption on 
the existence of moments. 

Definition 2.1. Let x be a random variable with (p{x) ^ 0. Then its 
S -transform is defined as follows. Let x denote the inverse under 
composition of the series 

oo 
n=l 

then ^ 

Sx{z) := x{z) • 

z 

The above definition has to be understood on the level of formal 
power series. Note that 'p)[x) ^ ^ ensures that the inverse of -0 exists 
as formal power series in z. Furthermore, can be recovered from Sx 
and thus the S'-transform contains all information about the moments 
of the considered random variable. 

The relevance of the S'-transform in free probability theory is due to 
the following theorem of Voiculescu [7] . 

Theorem 2.2. Ij x andy are free random variables such that ip{x) ^ 
and ip{y) ^ 0, then we have 

Note that in the above setting ip{xy) = ip{x)ip{y) ^ 0, thus the S'- 
transform of xy is also well-defined and the above theorem allows to 
get the moments of xy out of the moments of x and the moments of y. 
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Note that in this theorem we do not require x or ?/ to be real or 
positive random variables. In this formulation, the theorem is true in 
full generality, since it is essentially a statement aboTit mixed moments 
of free random variables. However, if one wants to extract a probabil- 
ity distribution out of the S'-transform, then one has to make special 
requirements concerning x and y. In particular, one would like to start 
with a situation where Lp{x"') are the moments of a probability mea- 
sure /X on M (i.e., x should be a real- valued random variable, i.e., a 
selfadjoint operator) and where are the moments of a probabil- 

ity measure on R (i.e., y should be real- valued random variable or 
selfadjoint operator). Since x and y do not commute, xy will not be a 
selfadjoint operator, thus it is not clear (and in general will not be the 
case) that (p{{xy)"') are the moments of a probability measure. How- 
ever, if X IS cl positive random variable, then ^/x makes sense and the 
operator \/xy^ is selfadjoint and has the same moments as xy. Thus 
in this case we know that the 9?((a;?/)'*) arc the moments of a probabil- 
ity measure on M; this is called the multiplicative free convolution of fi 
and V and denoted by /x IE z^. So the operation /x, i— > /x IE is defined 
for probabihty measures /x, i/ on R if at least one of them is supported 
on the positive real line. 

Often one restricts both factors to be supported in M+. In this case 
the above theorem for the S'-transform gives a solution to the problem 
of calculating \xMv. But there are also many interesting cases where 
only one factor is supported on R+. In this case the description of 
/X Kl z/ in terms of the S'-transform seems to have the following flaw. In 
order to be able to invert the power series '^(^z) = X^^i V'l^")-^" '^^^ 
needs a non-vanishing linear term, which means that the mean f{x) 
of the variable should be non-zero. For measures supported on R+ 
this is (with the exception of the uninteresting case x — 0) satisfied, 
but for measures supported on M the mean ^p{x) might be zero and in 
such a case it seems that the ^9-transform cannot be Tiscd to calculate 
/iKIz/. One should note that this is not an artificial situation, but covers 
such common cases like taking the multiplicative free convolution of a 
semicircle with a distribution supported on R+. One could of course 
try to approximate mean zero situations with non mean zero ones; 
however, the main purpose of our note is to point out that actually 
the failure of the S'-transform description in the mean zero case is just 
ostensivc; with the right interpretation, the above theorem can also be 
used in that case to determine fiM u. 

So let us consider the situation that yp{x) = for some selfadjoint 
element. Let us first see whether we still can make sense out of the 
relation Sxy{z) — Sx{z) ■ Sy{z) and whether this still allows to recover 
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all moments of xy uniquely out of the moments of x and of y. Again 
we exclude the uninteresting case that x — 0; hence we know that 
(p{x^) > (otherwise, by positivity, x would be equal to 0) and thus 
our series i^lz) starts with a multiple of z'^. This means that although it 
cannot be inverted by a power series in z it can be inverted by a power 
series in y/z. This inverse is not unique, but there are actually two 
choices (which correspond to choosing a branch of ^/z)^, however, this 
ambiguity does not affect the final result for the moments of xy, if y 
has non-zero mean. Let us be more specific about this in the following 
proposition. 

Proposition 2.3. (1) Let ip be a formal power series of the form 

oo 

i^iz) = ^ anZ^ = 0:22;^ + a^z^ H 

n=2 

with a2 ^ 0. Then there exist exactly two power series in \fz which 
satisfy 

i^ixiz)) = z. 

If we denote these solutions by 

X{Z) = Yl /^fc^'^' = + + /53V^' + • • • 

k=l 

and 

00 

x{z) = J2 ~^kz'/' = A + P2Z + h^z'' + • • • 
fc=i 

then their coefficients are related by 

(2) Let ip be as above and consider the two corresponding S-transforms 

S{z) = x{z) ■ ^ 
z 

and 

S{z) = x{z) ■ ^ 
z 

Then S{z) and S{z) are of the form 
and 
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where 

% = (-l)Sfe. 

(3) Let X be a random variable with vanishing mean and denote by 
and Sx the corresponding two S -transforms, as constructed in (1) 
and (2). Let y be another random variable, with non-vanishing mean 
so that its S-transform is of the form 



Sy{z) = 5kZ^ = 5o + 5iz + 52Z^ + 



fc=0 



Then the two series Sxy{z) :— Sx{z)Sy{z) and Sxy{z) :— Sx{z)Sy{z) are 
of the form 

wz ^ — ' 

^ fe=0 

and 



z 

k=0 



and related by 
If we put 

Xxy{z) = Sxy{z) ■ Y^jT^' Xxy{z) = Sxy{z) ■ ^ 

and denote by ipxy the unique solution of 

AyiXxyiz)) = Z 

and by ipxy the unique solution of 

'^'xyiXxyiz)) — Z, 

then we have 

i^xy(z) = -ii^xyiz). 

Proof. (1) The equation ■^(^(z)) = z means 



oo oo 



■+kn)/2 



n=2 k=l n=2 fei,...,fc„=l 

By equating the coefficients of powers of ^/z this is equivalent to the 
system of equations 

1 = a2Pl 

and 



= E E ■■■Pkn 



n=2 ki,...,kn=l 
fclH ^-kn=r 
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for all r > 2. 

The first equation has the two solutions 

ft ' 

and 



«2 



Pi = = -Pi 



"2 



Writing the other equations in the form 

r 

= 2a2Pr-iPi + Y1 Yl ^^"'^^i ■ ■ ■ -^^^ 

n=2 fei,...,fen=i 

fclH ^-kn=r 

(where the • • • indicate that we exclude in the second sum the cases 

ki = l,k2 = r — 1 and ki = r — l,k2 = 1) one sees that the values of /9„ 
for n > 1 are recursively determined by /9i and the a's. By induction, 
it follows that /3i = — /3i results in Pk = (— for all k. 

(2) This is clear, since we have 7^ = Pk+2+Pk for all k — —1, 0, 1, • • • , 
where we set := /?o := 0. 

(3) This follows by reverting the arguments from the first and second 
part. 

□ 

Definition 2.4. Let a; be a random variable with if>{x) = and <^(a;^) 7^ 
0. Then its two S-transforms S^; and are defined as follows. Let x 
and X denote the two inverses under composition of the series 

00 



n=l 



then 



Sx{z) := x{z) ■ and S^{z) := x{z) ■ 

z z 



Both Sx and Sx are formal series in y/z of the form 

yfc/2 



Now we stand ready to formulate our main theorem 

Theorem 2.5. Let x and y be free random variables such that (p{x) ~ 
0, f{x^) 7^ and ip{y) ^ 0. By Sx and Sx we denote the two S- 
transforms of x. Then 

Sxy{z) = Sx{z) ■ Sy{z) and Sxy{z) = Sx{z) ■ Sy{z) 
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are the two S -transforms of xy. 

Note that xy falls into the realm of Def. I2.4[ We have 

(p{xy) = (p{x)ip{y) = 

and 

ipiixyY) = ifix^Myf + v'(x) V(?/') - vix)\{yf = v^i^^Myf 7^ 0. 

Let us point out again that the ambiguity of having two ^'-transforms 
for xy is not a problem, because both of them will lead to the same 
moments for xy, according to Prop. 12.31 

In order to prove our Theorem 12.51 one might first inspect the usual 
proofs for the S'-transfom. One notices that all of them rely on the fact 
that (p{x) and ip{y) both are not zero. None of the published proofs can 
be used directly for the case where one variable is centered. However, 
at least the idea of the proof in [3l H] can be adapted to our situation. 
The following proof is a variant of the approach in [Sj Hj (by avoiding 
an explicit use of the incomplete boxed convolution) and can of course 
also be used to give a straightforward proof in the usual case (p{x) ^ 0. 

Proof. For a random variable x we denote by 

oo 
n=l 

and 

M,(z) := l + 
the corresponding moments series and by 

oo 
n=l 

the corresponding cumulants series (k^ = Kn{x, . . . ,x) is here the n-th 
free cumulant of x). 

Consider now x and y as in the theorem. Then, in addition to the 
moment and cumulant series of x, y, and xy, we consider also additional 
moments series of the form 

oo 
n=0 

and 

oo 

M,{z) :=Y,if{x{yxr)z\ 

n=0 

By the moment cumulant formula and by the fact that because of the 
freeness between x and y mixed cumulants in x and y vanish (see [1] 
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for the combinatorial theory of freeness), it is quite straightforward to 
derive the following relations between these power series: 

(1) M^y{z) = My,{z) = Cy[zM2{z)] + 1 

M,y{z) 



(2) Mi(z) = Cy[zM2{z)] 

(3) M2{z) = C^[zM^{z)]- 



zM2{z) 
M^y{z) 



zMi{z) 

Note that all these relations are valid (and make sense as formal power 
series) independent of whether if{x) = or not. 

Let us note that we have the well-known relation M{z) = 1 + 
C[2;M(2;)] between moment and cumulant series, which shows, by re- 
placing z with x{^)y that C[^S'(2;)] = z. This relation, which was 
observed in 0, can be taken as an alternate definition of S{z). In the 
case where (p{y) 7^ this is just a relation between formal power series 
in z. In the case f{x) = 0, it has again to be read as formal power 
series in ^/z. (Note that for the case ip{x) = 0, (/?(x^) 7^ we also have 
= if{x) = and Kg — V^l^^) 7^ 0; thus Cx{z) starts also with a 
quadratic term in z.) 

Since f{xy) = (f{x)ip{y) = 0, the moment series ipxy has two inverses; 
let X be one of them. Replacing z by x{^) dl]) gives 

l + z = M^yixiz)) = Cy[x{z)M2{x{z))] + 1, 

yielding that 

(4) z = Cy[x{z)M,{x{z))]. 
Equation ([2]) gives 

zM^{z)M2{z) = Cy[zM2{z)] ■ M^y{z). 
Replacing z by x{^) gives 

X{z) ■ M^{xiz)) ■ M^ixiz)) = Cy[xiz)M2{x{z))] ■ M.,(x(^)) 
and hence 

(5) x{z) ■ Ml {x{z)) ■M2{x{z))=z-{z + l) 

Since zSy{z) is the unique inverse of Cy under composition we must 
have, by that 

zSy{z) = x{z)M,{x{z)). 
In the same way as for (jl]) we get 

z = C4xiz)M,ix{z))]. 
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Since Cx has exactly two power series in ^/z as inverses under compo- 
sition, namely zSx{z) and zSx{z) we must have either 

zSx{z) = x(z)Mi(x(^)) or z~Sx{z) = x(^)Mi(x(^)). 

Name 5*2; and Sx in such a way that we have the first equation. 
Plugging this into ^ gives 

zSx{z) ■ zSy{z) = xiz)MMz)) ■ x(^)M2(xW) = xiz)ziz + 1) 

and thus 

z + 1 

Sx{z)Sy{z) = X{.Z)^Y~ = ^xy{z). 

If we start with the other inverse, x? of i^xy, then we would end with 
the other S'-transform, Sxy □ 

This proposition tells us that we can also use the formula Sxy{z) = 
Sx{z)Sy{z) to calculate the moment series of xy, for x and y free, in the 
case where at most one of them has non-vanishing mean. Since this 
is the case in the situation where y has a distribution supported on 
M^, we see that the S'-transform is a useful tool for the calculation of 
/i Kl zy, whenever /i is a probability measure on M and u is a. probability 
measure on M_(_, independent of whether /i has vanishing mean or not. 

3. Case where both means are zero 

One should also note that for the proof of our proposition it is impor- 
tant that at least one of the two involved measures has non-vanishing 
mean. If both have vanishing mean then the arguments break down. 
One should, however, note that this cannot be attributed to the fact 
that there does not exist a probability measure with the moments of 
xy. Actually, for x,y free with ip{x) = = (p{y) we have by the def- 
inition of freeness that (p{{xy)"') = for all n, i.e., xy has the same 
moments as the delta distribution 60- So one might say that the mul- 
tiplicative free convolution of any two measures with mean zero is Sq. 
However, on a formal level the situation for the S'-transform gets quite 
different here. As a concrete example, let us consider the multiplicative 
free convolution of the semicircle with some u. The S'-transform of the 
semicircle (of variance 1) is 

S^,{z) = 



If we multiply this with an S'-transform which is a power series in z 
then the result is again a series of the form 

■ power series in V^, 



10 



N. RAJ RAO AND R. SPEICHER 



which can, as in our proposition, be uniquely resolved for the corre- 
sponding moment series. Let us take now, on the other hand, the 
multiplicative convolution of two semicirculars. Then we have 

Sf,{z)Sf,{z) = - 
z 

which results in a transform of the form 

,p{z) = ^ - 1. 

This, however, is not a valid moment series. On the other hand, 60 
would have the moment series i^iz) = 0, for which no corresponding S- 
transform exists. So we see that in the case where we multiply two free 
operators with mean zero the S'-transform machinery breaks totally 
down. 



4. Some examples 

Consider the free multiplicative convolution of the semi-circle distri- 
bution with the free Poisson distribution. The semi-circle distribution 
has S'-transform 



'z 

while the free Poisson distribution has S'-transform 

S^(z^ 



z+l 

The distribution obtained by their free multiplicative convolution has 
S'-transform 

^'^^^^^^ = ^z{z + iy 

The Cauchy transform of the probability measure Kl 7 satisfies the 
algebraic equation 

g^z^-zg + 1 = 0, 



from which we can obtain the density function, plotted in Figure l(a' 



Consider now the free multiplicative convolution of the free Poisson 
distribution with the shifted (by —a) free Poisson distribution. The 
free Poisson distribution (shifted by —a) has S-transform 



, -z-l + a + j^2z + 2za + 1- 20 + 0"^ 

^^^'^ = 

which, for a = 1 yields 



^ , . ^^ , . -Z + -J Z^ +4:Z 
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The distribution obtained by the free multiphcative convolution of the 
free Poisson distribution with the free Poisson distribution (shifted by 
— 1) has (S-transform 

^''^^'^ - 2z(l + .) • 

The Cauchy transform of the probabihty measure yU Kl 7 satisfies the 
algebraic equation 

g^z'^ + z'^g^ - zg"^ - gz+1 = 



from which we can obtain the density function, plotted in Figure 1(b) 
We note that the computations in this section were done using RMTool 
software |5j based on the algebraic equation based framework for com- 
putational free probability developed in [6l [2] . 
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(a) The solid line is the density function of the probability measure obtained by the 
free multiplicative convolution of the semicircle with the free Poisson distribution. 
The histogram bars are the eigenvalues of the product of a Wigncr matrix with a 
Wishart matrix and were generated using 50 x 50 sized random matrices over 4000 
Monte-Carlo trials. 
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(b) The solid lino is the density function of the probability measure obtained by 
the free multiplicative convolution of the free Poisson distribution with the shifted 
(by —1) free Poisson distribution. The histogram bars are the eigenvalues of the 
product of a Wishart matrix with an independent Wishart matrix (shifted by the 
negative identity) and were generated using 50 x 50 sized random matrices over 4000 
Monte-Carlo trials. 
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